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1. Introduction 



Recent astrophysical observations of distant galaxy distributions |T]|,[|1 and the cosmic 
microwave background have revealed a quite astonishing picture of the universe. From the 
survey of relatively near galaxies, the void structure (the Great Wall) was discovered, p| and 
from the pencil-beam survey of galaxies, the quasi-periodic distributions (of period about 
ISOMpc) were inferred. [|l| The data from COBE, on the other hand, have revealed extremely 
isotropic and homogeneous distribution of the 2.7 K cosmic microwave background with fluc- 
tuations of order 1O~^.0] In general, it is very difficult to explain how these anistopic and 
inhomogeneous large-scale structures of the universe have developed from such an isotropic 
and homogeneous distribution of matter in the early stage of the universe. The standard 
solution to this difficulty totally relies on the existence of dark matter which accounts for 
more than 90% of matter in the universe. The evidence for dark matter was first claimed in 
order to explain the flat rotation curves of spiral galaxies. Since there is no established direct 
observation of dark matter, however, there are many attempts to explain the rotation curves 
without dark matter by modifying the Newtonian force 0] or by modifying the Newton's 
second force law.||^,|]^ Other people have tried to derive such modifled Newtonian force laws 
from the framework of general relativity. [|] [1^ 



In a previous paper, |T0| we attempted to explain not only the flat rotation curves of spiral 
galaxies but also the large-scale structure of the universe, starting from a simple model with 
the addition of a quadratic scalar curvature term to the Einstein action. Our generalized 
action could qualitatively explain the flatness of the rotation curves and the nearly periodic 
galaxy distributions. However, it turned out that our theory does not imply the unity of the 



coefficient 7 of the Robertson expansion [|TT| on the solar scale. This constraint (7 = 1 



from the classical tests of general relativity such as the observation of the radar echo delay 
is quite stringent, and it is very difficult to realize this value in the generalized Einstein action. 



In this paper, we construct the generalized Einstein action under the two criteria: i) it 
must give 7 = 1 in the post Newtonian approximation, ii) the gravitational potential is a 
sum of Newtonian and Yukawa potentials. The second criterion is imposed, since it is the 
empirical gravitational potential of Sanders that can quite successfully explain the flat 
rotation curves of spiral galaxies. We then show that the minimum ingredient of the theory 
that satisfles the above criteria is the term, a scalar fleld, and a vector fleld in addition 
to the Einstein action. 
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2. Generalized Einstein Action and Its Post Newtonian Approxi- 
mation 



The generalized Einstein action which contains quadratic terms of the scalar curvature, 
i?^, and Ricci tensor, R^^R^" , was introduced to regulate the ultraviolet divergences of the 
Einstein theory. [p!2| It was applied to cosmology to obtain the bounce universe to avoid the 
singularity at the creation of the universe. [13| The structure and the properties of the theory 
were further elaborated in subsequent works. WA 



In this section, we consider a further generalization of the theory by adding scalar and 
vector fields in addition to R^ and R^^R'^'^ terms and study its post Newtonian approxima- 
tion. We will investigate such a theory with two requirements: i) on the solar scale, it must 
conform with the classical tests of general relativity, ii) on the galactic scale, the gravita- 
tional potential is modified to give a sum of Newtonian and Yukawa potentials in order to 
explain the rotational velocity curves of spiral galaxies. 

We consider the following generalized action with the scalar curvature R, the Ricci tensor 
R^^, a scalar field a, and a vector field A^: 

I = J d'xV^{-Y^{R + c,R' + c^R^^R,,) 

-U,ad^a-^a' + g,aR (1) 



where G is the gravitational constant, (i and m are the masses of the scalar and vector par- 
ticles, is the covariant derivative, and is the matter Lagrangian. The coefficients 
ci,C2, and G have the dimension of (mass)"^, while the coefficient gi has the dimension of 
mass and (72 is dimensionless. 



In order to calculate the coefficient 7 in the Robertson expansion, ||11|1 we introduce the 
weak fields and ifj defined by 

(700 = -1 - 20, g^j = 6ij{l + 2ij) . (2) 

In the post Newtonian approximation, we must take into account up to the quadratic term 
of the weak fields and source in the action, and the necessary formulae to the first order in 
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the weak fields are 



Rij ~ didj{(j) + ^) + 6,jA^ , i? ^ 2A(0 + 2iIj) 



(3) 



and the necessary formulae up to the second order in the weak fields and source in the action 
are 

^/^R ^ 2A(0 + 2^) + 4^ Alp + 2ipAip , 



^R^'^R^^ ^ 2(A0)2 + 4A0A'(/' + 6{A^f 



(4) 



-gR^^A{A{<j) + 2^)y , 



matter 



where we have suppressed the total derivative terms. ||15| We substitute the weak field ex- 
pressions Eqs. (^) and (|p into Eq. (|T]), and retain up to the quadratic terms of the weak 
fields 0, if), 0", and source p to obtain 



d X 



1 



16nG 



{2A(0 + 2ip) + A(f)Aip + 2ipAip + 4ci(A(0 + 2ip)y 
+4c2(i(A0)2 + A0Az^ + ^(A^)2} 



1 

— L 

2 



+-aAa - ^a' + ^AAA, - '-^A^ - ^AqAAo + '-^Al 



m 

T 



m 

Y 



(5) 



+2giaA{(P + 2^) + 2^25iAA(0 + 2^^) - p0 

Since the field Aq decouples from the other fields and source, we do not consider this field 
hereafter. It is convenient to introduce the following variables and ip: 



which in turn gives 



+ 2^ , ^ = ^-^ 



- 27/^ , + V 

— ^ — > = 



(6) 
(7) 



Substituting Eq. (|^ into Eq. we can rewrite the action as 



1 



,_1_{2A0 + ^0A0 - + 4ci(A0)2 + ^((A0)2 + i(A^)2)} 



P0 , 2p^ 



-aAa - —a^ + -AiAAi Af + 2^iaA0 + 2g2diA^A(j) - ^ + 



(8) 



Taking the variation with respect to the weak fields 0, ip, Ai and a, we obtain the equation 
of motion of the weak fields in the post Newtonian approximation: 



A {1 + (6ci + 2c2) A} (p - 2ATxG{giAa + g2Ad,Ai) = -AuGp 



(9) 



A(1-C2A)V' = -SvrGp 



(A - m2)Ai - 2g2Ad,(f) = 



(10) 



[A - fi^)a + 2giA(t> = 



:i2) 



From Eqs. ([TT|) and ([l^), we have 

2^2AV 2^iA0 
A — A — fi^ 

Substituting this expression into Eqs.(P) and (p!OD , we then obtain the following: 

A (l + (6ci + 2c2)A - A8nG{--f^ + -^^A = -^nCp , (14) 
[ A — A — J 

A{1- C2A)ip = -8nGp . (15) 



In the next section, we will consider the possibility that the coefficient of the Robertson 
expansion 7 becomes 1. In order to obtain this result, it is necessary that both (p and ip 
behave ~ 1/r in the limit r —>■ 0. { We remark from Eq. (|^): if and ip> behave like (p ~ const, 
and ip 1/r in the limit r — > 0, we obtain 7 = 1/2, while if (p and ip behave like ~ 1/r 
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and ~ const, in the limit r — 0, we obtain 7 = —1.) In the region r — 0, the mass term 
is neghgible, and the necessary condition to have 7 = 1 is that the highest derivative terms 
(cx A^) on the left-hand sides of Eqs. (|Iip and ( [T3| ) must vanish (for fi^ = m? = 0). This 
condition reads 

6ci + 2c2 - ASirGgl = , Ca = . (16) 
Near the origin we then have 

(1 - Gcim^ + A8nGgl)A(p ^ -AnGp , (17) 

Aij ^ -8nGp . (18) 



We assume that the density takes the point-like distribution of the form 

p(f) = M6{r} . 
Using the formula AnGp/A = —GM/r, we then obtain 

- _ AnGp _ GM 
kA kr ' 



(p-2'ip _ GM Ak-1 GM 
3 r 3k r 



(19) 



8nGp 2GM , , 

^^-^ = --, (20) 



where 

k = l- 6cim^ + ASTiGgl . (21) 
Therefore we obtain the gravitational potential of the form 



f22) 



<P + ^P GM 2k + l GM 



where 
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This leads to the following formula for 7 : 

7 = ^. (25) 

Here, by taking Morikawa model as an example, we demonstrate how it is difficult to 
satisfy the stringent condition 7 = 1 in a modified Einstein theory in general. 

Morikawa model is a Brans-Dicke type theory of the form |[T^ 

-\d,'fd^'f - f^v' + >-V + + i™""} , (26) 



rMorikawa 



where we have added the cosmological term A and tadpole term Xcp to the original Morikawa 
action, since we consider the case that the scalar field ip takes an expectation value v. We 
write (f = V + a and consider this a field as the weak field. We then obtain the weak field 
approximation of the form 



rMorikawa _ [d^^^/^I —h^SnG^V^' 

J ^ 16-kG 



1 

--d^ad^'a - + ^vaR + L™<^«-"| . (27) 



Here, we have tuned A and A in such a way that 

A li^v'^ 



+ Xv) = , A - /i^t; = 



for a given v. If we denote 1/G = (1 — 8'n-G^v'^)/G, then according to our formula Eq. pT| ) 
(with the replacements m — >• 0, G ^ G, and gi ^v), we have the expression k = 1 + 
ASirG^'^v'^, which gives 

l-8.G(e-4eV_ (28) 

In Morikawa's case, [0 the scale factor of the universe oscillate with time, where the period 
is converted to the scale of ISOMpc through the velocity of light, and it is unclear whether ip 
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here takes an expectation value on the solar scale. We obtain 7 = 0.88 by using Morikawa's 
unit 1/G = Air/S and his values ^ = 10 and v = 0.008, where the value of v is assumed to 
be of the order of the initial value of (f . (Of course, if we set f = in Morikawa model, we 
have 7 = 1.) 

Through this example, we understand that the construction of the generalized Einstein 
action with 7 = 1 is quite non-trivial. 



3. Generalized Einstein Action on Solar and Galactic Scales 

In this section, we show that the generalized Einstein theory, in which the gravitational 
potential is a sum of Newtonian and Yukawa potentials, must include at least scalar and 
vector fields in addition to the quadratic term of the scalar curvature in order that the con- 
dition 7 = 1 is satisfied. 

The condition 7 = 1 implies = 1 in Eq. (^), which in turn implies Cim^ = SnGgf 
(see Eq. (pT])). Taking into account Eq. (|T6|) also, the necessary condition to have 7 = 1 is 
summarized as follows: 

ci = SttCqI, C2 = , (29) 
cim^ = STTGgl . (30) 



This condition is classified into the following possibilities: 

/) 92 7^0, 

a) gi = , (ci 7^ 0, m = 0), 

b) gi^O, (ci ^ 0, m ^ 0), 
//) 92 = 0, (ci = 0, gi = 0). 

We first consider case la) , where the scalar decouples from other fields. In this case Eqs. ([I^ ) 
and ( jT3| ) become 

A<^ = -AnGp, Atp = -8nGp , (31) 

which in turn gives 

AnGp GM - SnGp 2GM , , 

^ = ^ = — • (32) 



A r ' ^ A 
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Hence, we obtain the same result as the Newtonian approximation of the ordinary Einstein 
theory: 

^ = — ' ^ = ' (33) 

in which the desired Yukawa term is absent. 

In case IT), the vector field decouples and there is no higher derivative term such as 
and Rfj,uR^^ . In this case, Eq. (|1^) takes the same form as Eq. (0), and we obtain the result 
of the Newtonian approximation again. 

The final case Ih) turns out the one that satisfies our criteria. In this case, from the 
relation cim? = SnGgf, Eq. (|T^) becomes 

___ji_(A'^ - ^ ('"^ + + - ,^)) 

+mV}0 = -AnGp . (34) 

We then define a and P by 

a + f3 = + + AMGglim^ - p^) , 

a(3 = m^/i^ , (35) 

and assume a > (3. We can also define constants ki, k2 and by 

- 47rG'(A-m2)(A-/i2) 

-P , 



A(A-a)(A-/3) 



which gives the relations among ki, k2 and k^ as follows: 

^1 + ^2 + ^3 = 1 , 

ki{a + (3) + k2l3 + k^a = + /x^ , (37) 
kial3 = m^i/ . 

From Eq. (^) we obtain fci = 1, so that k2 and k^ are related by k^ = —k2- Using this 
relation and substituting Eq. (^) into Eq. (P7|), we obtain 

487rG^f (m^ - /i^) 

«;2 = -ks = ^ • (38) 

a — p 
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The modified gravitational potential in this case becomes 

GM 



1 + k2e-^'' - k2e-V^A , (39) 



r 

2GM 

ijj = . 

r 

Therefore we have the gravitational potential of the desired form: 



We now examine the condition for this gravitational potential to explain the flat rotation 
curves of spiral galaxies. We know that the gravitational potential in Sanders' form [ffl 



GM fl + as e~ 



1 + ^5 



(41) 



can account for the rotation curves in a satisfactory way when as = —0.9. For our potential 
Eq. (^OD to have a similar form on the galactic scale, it is necessary to assume (for a > f3 > 0) 

a > f3 > , (43) 
where Vq is a distance on galactic scale (~ a few lOkpc). Note that these relations imply 

exp(-y^ro) ^ 1 . (44) 

For the rotation curves to be flatter than the Newtonian results, we need < in Eq. (^Tj). 
In our potential Eq. (HO) this requirement together with Eqs. (H3|) and (HI) implies 



^2 > . (45) 

By Eq. (|38|) this is equivalent to 

> /i^ . (46) 
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Hence, the question is whether we can satisfy the above conditions Eqs. P^), (^31), and 
(or (46)) in a consistent way. From Eqs. (|35|) and (^61), one obvious solution is 



= 0{f,') , (47) 

where ro is the galactic scale. In this case, from Eq. (|35D a can be written as 

a ^ m^{l + 487iGgl) . (48) 
By defining x = 4:87!'Ggf{> 0) and neglecting /? and fj,"^ in Eq. (pS]), we have 

so that /c2 can change in the range 

< A;2 < 1 . (50) 
The coefficient of Eq. (|4l|) in our theory is then written as 

h X 

6 + k2 6 + Ax 

Note that our as takes a value in the range 

- 0.25 < 05 < , (52) 
while Sanders takes as = —0.9 to fit the rotation curves. 

In this way, starting from the generalized Einstein action with additional scalar and 
vector fields, we can tune the parameters so that we not only have 7 = 1 on the solar scale 
but also obtain the empirical Sanders type gravitational potential on the galactic scale. 



4. Summary and Discussion 

In this paper, we have examined the generalized Einstein theory which contains higher 
derivative terms, and R^'^R^^, and satisfies the criteria: i) on the solar scale, it must be 



11 



consistent with the classical tests of general relativity, ii) on the galactic scale, the gravita- 
tional potential is a sum of Newtonian and Yukawa potentials so that it may explain the flat 
rotation curves of spiral galaxies. We have shown that it is non-trivial to satisfy the above 
criteria and that at least additional scalar and vector fields are required for a consistent 
theory. 

We have tuned the parameters of the theory so that the coefficient as of the Yukawa term 
( in Sanders' gravitational potential) is negative, which is necessary to explain the flatness of 
rotation curves of spiral galaxies. It will be interesting to see how well our theory can fit the 
rotation curves quantitatively. Numerical calculations in this direction are now in progress. 

In our generalized Einstein action, even after the tuning of parameters, there are still two 
parameters, the vector mass m and the scalar /x, to set the scale of interest. We have chosen 
1/m ~ galaxy scale, and it may be possible to explain the "periodic" large-scale structure 
of the universe by choosing ~ ISOMpc. Numerical calculations of distant galaxy distri- 
butions are also in progress. 

Therefore, by taking the values of 1/m and 1/n as above, our generalized Einstein theory 
may be consistent with observations over three different distance scales: the solar, galactic 
and beyond galactic scales of the universe. 
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